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Abstract 

We review some aspects of the isotope effect (IE) in superconductors. Our 
focus is on the influence of factors not related to the pairing mechanism. After 
summarizing the main results obtained for conventional superconductors, we re- 
view the effect of magnetic impurities, the proximity effect and non-adiabaticity 
on the value of the isotope coefficient (IC). We discuss the isotope effect of T c 
and of the penetration depth 5. The theory is applied to conventional and high- 
T c superconductors. Experimental results obtained for YP^CuaOy^ related 
materials (Zn and Pr-substituted as well as oxygen-depleted systems) and for 
La 2 _ x Sr a; Cu04 are discussed. 



1 INTRODUCTION 

Historically, the isotope effect (IE) played a major role in unravelling the questions 
related to the origin of the effective attractive interaction between charge-carriers which 
leads to the superconducting state. The theoretical considerations of Fr6hlich|l[ and the 
experimental discovery of the IE of T c in mercury^ pointed towards the contribution 
of lattice dynamics to the instability of the normal state. Since then, the IE has often 
been considered as a measure of the contribution of phonons to the pairing mechanism. 
Furthermore, it is generally assumed that only those thermodynamical quantities that 
depend explicitely on the phonon frequencies (as T c or the order parameter A in the 
BCS model) display an IE. 

One can show|3|, §, || @, |7j, however, that such understanding of the isotope effect 
is incomplete and leads, therefore, to confusion; this is particularly true for the high- 
temperature oxides. Indeed, several factors not related to the pairing mechanism can 
alter the value of the isotope coefficient (IC). Moreover, these factors are not necessarily 
related to lattice dynamics. Here we focus our attention on three such factors: magnetic 



impurities, the proximity effect and non-adiabatic charge-transfer as it occurs in high-T c 
superconductors. 

In addition, we also show || that fundamental quantities such as the penetration 
depth 5 of a magnetic field also display an IE because of the three factors mentioned 
above. Note that this effect occurs despite the fact that 5 does not explicitely depend 
on quantities related to lattice dynamics. 

It results from our considerations that the value of the IC (even its absence) does not 
allow any a priori conclusion about the pairing mechanism. Nevertheless, it remains 
an interesting effect that enables one to determine the presence of magnetic impurities, 
proximity effect or non-adiabaticity in the system. The calculation of the IC and its 
comparison with experimental results was also used in previous works to determine the 
value of the Coulomb repulsion /z* or the relative weight of different electron-phonon 
coupling strengths in superconductors (see below). It can thus be used as a tool for 
the characterization of superconductors. 

This review is mainly based on our papers 0, ||, [5], §, We apply the theory to 
analyse the oxygen isotope effect of T c in Zn and Pr-doped (YBCZnO and YPrBCO) as 
well as oxygen-depleted YBa 2 Cu307_,5 (YBCO). We also review calculations of the pen- 
etration depth isotope effect. The latter superconducting property is a good example 
of a quantity that does not directly depend on phonon frequencies and yet can display 
a substantial isotopic shift ||. In this context, we also discuss recent experimental 
results obtained for La2- x Sr x Cu04 (LSCO). 

The structure of the paper is as follows. In section ^| we present some of the early 
results related to the isotope effect and show that even for conventional superconductors 
the relation between the isotope effect and the pairing mechanism is not simple. We also 
discuss the applicability of these early results to the description of high-temperature 
superconductors. 

The remaining sections are devoted to the description of new, unconventional isotope 
effects. In section [3] we study the influence of magnetic impurities on the value of the 
isotope effect. We show that adding magnetic impurities to a superconductor can 
enhance the isotope coefficient a of T c and induce a temperature-dependent isotope 
effect of 5. We show that both isotope effects are universal functions of T c . That is, 
a(T c ) and (3{T C ) are independent of any adjustable parameter. We discuss Zn-doped 
YBa 2 Cu 3 07_,5 in this context. 

Section ^ is concerned with the influence of a normal layer on the isotope effect of a 
superconductor. We show that due to the proximity effect, the isotope coefficient of T c is 
linear in the ratio of the normal to the superconducting film thicknesses p. Furthermore, 
the proximity effect induces a temperature and p-dependent isotope coefficient of the 
penetration depth. 

Section |^ reviews the concept of the non-adiabatic isotope effect introduced in 
Ref. and further discussed in Refs. (|, |5|, ||, We show that in systems as the 
high-temperature superconductors where charge-transfer processes between reservoir 
and Cu02-planes occur via ions that display a non-adiabatic behaviour, the charge- 
carrier density in the planes depend on the ionic mass. This leads to the unconventional 
non-adiabatic contribution to the IE of T c . It is also interesting that non-adiabaticity 
induces an isotopic shift of the penetration depth S. 

In the last section (Sec. H) we apply our theory to the oxygen isotope effect (OIE) 
in high-temperature superconductors of the YBCO-family. Focus is set on the oxygen 
isotope effect because most experimental studies have been performed on the oxygen 
ion which is the lightest in the Cu02 plane. The situation is less clear in the case of 
copper isotopic substitution where an effect has also been observed []8], [|. This latter 



case will be discussed in more detail elswhere. 
We conclude the review in section [?]. 

2 THE ISOTOPE EFFECT: DISCOVERY, 
CONVENTIONAL VIEW 

The isotope effect of superconducting critical temperature T c is best described in terms 
of the isotope coefficient (IC) a defined by the relation T c ~ M~ a , where M is the 
ionic mass. Under the assumption that the shift AT C induced by isotopic substitution 
(M — > M*) is small compared to T c , one can write 

M AT C 



where AM = M — M* is the difference between the two isotopic mass. If the supercon- 
ductor is composed of different elements, one defines a partial isotope coefficient a r as 
in Eq. ([!]), but where M is replaced by M r , the mass of element r that is substituted 
for its isotope. In addition, one defines the total isotope coefficient by 

Besides T c there are other quantities that display an isotope effect. In the next 
sections we focus on the isotope shift of the penetration depth 5. In analogy to the IC 
of T c (denoted a) we define the isotope coefficient (3 of the penetration depth by the 
relation 

M A5 

One should note that there is a conventional BCS-type isotope effect of the penetration 
depth related to its temperature dependence 5~ 2 ~ (1 — T 4 /T c 4 ). Because T c displays 
an isotope effect, the penetration depth is also shifted upon isotopic substitution. The 
effect becomes strong as one approaches T c . The present article is not concerned with 
this trivial effect. 

Let us first discuss what values of the isotope coefficients a and (3 can be expected 
for different types of superconductors. Table 1 shows characteristic values of the isotope 
coefficient for different types of superconductors. Very different values of the IC have 
been observed (in the range from —2 to +1). One notes that some systems have a neg- 
ligible coefficient, some display even an "inverse isotope coefficient" (a < 0) and some 
take values greater than 0.5, the value predicted by Frohlich and the BCS model (for a 
monoatomic system). The purpose of the next paragraphs is to describe shortly differ- 
ent theoretical models allowing one to understand the coefficients observed (see Table 1 



and Refs. [ IP], [TT|] ) . We begin with the description of conventional superconductors and 
discuss the relevance of these models for high-T c materials. In the following sections 
we then introduce new considerations about the IE allowing one to give a consistent 
picture of the IE in high-temperature superconductors. 



Table 1. Experimental values of the isotope coefficient of T c (see also Refs. [IT , |20"fl). 
The letters in the last column correspond to Ref. || . 
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2.1 Monoatomic Systems 

Let us begin with the case of a monoatomic BCS-type superconductor. Substitut- 
ing the atoms by an isotope affects the phonon dispersion (for a monoatomic lattice 
f2 ~ M -1 / 2 , where £1 is a characteristic phonon frequency). Thus, any quantity that 
depends on phonon frequencies is affected by isotopic substitution. Assuming that 
the electron-electron pairing interaction is mediated by phonons and neglecting the 
Coulomb repulsion between electrons, the BCS theory predicts that T c oc fl (see below, 
Eq. (^) for fi* = 0). Thus, for a monoatomic system the isotope coefficient of T c given 
by Eq. (0) is a = 0.5. 

As seen in Table 1, the IC of most monoatomic non-transition metals is approxi- 
mately equal to 0.5. Many other systems, on the other hand, deviate from this value. In 
particular, transition metals and alloys display values that are smaller than, or equal to 
0.5. PdH displays an inverse isotope effect ctpdu — —0.25 |H| (there is also one report 
on a large inverse IC of Uranium otjj ~ — 2|12|]). Furthermore, high-T c materials can 
have values both smaller and larger than 0.5, depending on the doping. 

2.2 The Coulomb Interaction 

One of the first reasons advanced to explain the discrepancy between theory (a = 0.5) 
and experiment was that the BCS calculation described above did not take properly 
into account the Coulomb repulsion between charge-carriers. It was shown in Refs. |13| , 
[TR [TB| that inclusion of Coulomb interactions leads to the introduction of the 
pseudo-potential /i* in the BCS equation for T c : 

T c = 1.13ft exp^-^-^) , (4) 

where /i* is the pseudo-potential given by 



/i is the Coulomb potential and Ep is the Fermi energy. One notes that ft is present 
in /i* and thus in the exponent of Eq. (Q). The value of a can thus be substantially 
decreased from 0.5. Indeed, from Eq. (j^) and (f|) one obtains 



a 



(6) 



The isotope coefficient is shown as a function of A for different values of /i* in Fig. 1. 
Naturally, only the case A > //* is relevant, since superconductivity is otherwise sup- 
pressed. 

The effect of Coulomb interactions can be included in a similar way in the strong- 



coupling Eliashberg theory. Using the formula derived by McMillan [17 




(7) 



with fj, = (1 + 0.62A)/i* and /i* is given by Eq. (H), one obtains the following result for 
the isotope coefficient: 
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(8) 



A similar result has also been derived in Ref. Fig. 1 shows the dependence described 
by Eq. (^) and allows one to compare the weak and strong-coupling cases, that is, 
Eqs. (§) and (g). 

An expression for T c that is valid over the whole range of coupling strengths (from 
weak to very strong couplings) has been derived by one of the authors in Ref. |H 



0.25ft 



with 



yJ&TKff - 1 

1 + 2/x* + Xfj*t{X) 

and t(X) ~ 1.5exp(— 0.28A) |[T8|] . As usual, the parameter /x* is given by Eq. 
isotope coefficient resulting from Eq. ([!]), (|^) and ( |I0"D reads 
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This isotope coefficient is also displayed in Fig. 1 for two values of the parameter /j*. 
Note that for all models, because of the presence of //* the isotope coefficient of T c 
is lowered with respect to the BCS value for a = 0.5. One recovers the BCS result 



asymptotically at large A (or T c /Vt). Only Eq. fllT]) , however, is valid in the strong 
coupling limit (A > 1). 

An interesting general feature of the results presented above is that the strongest 
deviation from the BCS monoatomic value 0.5 occurs when /x* (or jX) is of the order 
of A (see Eqs. (^), (|8|) or ([□])). Thus, excluding anharmonic or band-structure effects, 
a small isotope coefficient is correlated to a low T c . This is indeed observed in most 
conventional superconductors. In certain cases one even obtains a negative (also called 
inverse) isotope coefficient (a < 0). For realistic values of the parameters, T c should not 




T c /n 



0.15 



Figure 1. Isotope coefficient a as a function of T c /Q. (i.e. A) for different values of /x*. Solid 
lines from Eq. (||), dashed lines from Eq. (^) and dash-dotted lines from Eq. (11). For each 
type of line, the upper curve corresponds to u* = 0.1 and the lower to //* = 0.2. 



exceed T c ~ IK for this inverse effect to be present. Several systems display this inverse 
isotope effect, among them PdH (where H is replaced by its isotope D) [[!3|, Q and 
Uranium (where 235 U is replaced by 238 U) [12] ■ We discuss this case below (sec. |2.3| ). 



The situation encountered in high-temperature superconductors is very unusual. 
Indeed, one observes that the isotope coefficient has a minimum value at optimal doping 
(highest T c ). The explanation of a small a for oxygen substitution coinciding with a 
high T c is likely to be related to the contribution of the oxygen modes to the pairing. 

The main conclusion of this section is that the Coulomb interaction and its logarith- 
mic weakening leads to the deviation of the isotope coefficient from the value a = 0.5 
and to the non- universality of a. 

2.3 Band Structure Effects; Transition Metals 

Up to now we have only considered the effect on T c induced by the isotopic shift of 
phonon frequencies. Furthermore, we have considered the case where the phonon fre- 
quency appears explicitely in the expression for T c (as prefactor and in /i*; see Eqs. (^), 
(0) and ©). Band structure effects have first been considered to explain the discrepan- 
cies observed between the results following from the strong-coupling McMillan equation 
(Eq. (0)) and the IC measured in transition metals. Several other generalizations of 
the basic model have been considered, especially in the context of high-temperature 
superconductors. 

Let us first describe the situation encountered in transition metals. As has been 
shown in Refs. |15|, pl| , the two-square model for which (|5|) is derived leads to values of 
the IC that are only 10 — 30% lower than 0.5. However, the deviation is much stronger 
for some transition metals, leading even to a vanishing coefficient for Ruthenium. As 
shown in Ref. [113 it is necessary to take into account the band structure of transition 



metals. The presence of a "metallic" s-band and a narrow <i-band, and the associated 
peaked structure of the electronic density of states (DOS) has an impact on the value of 
/i*. Furthermore, Ep is much smaller in transition metals than in non-transition metals 
and leads thus to a smaller effective screening of the Coulomb interaction. These facts 



have been included in Ref. [19j by modifying the two square-well model used to derive 
Eq. (H). The values obtained from the modified expression of a (see Eq. (9) in Ref. |19|| ) 
are in good semi-quantitative agreement with the experimental results. 

One should add that the two-band structure is important even in dirty transition- 
metal superconductors, although the superconducting state is characterized by only 
one gap in this case (see Ref. |L1| for further discussions on transition metals). 

The example of transition metals shows that a quantitative description of the IC 
requires, among other things, a precise knowledge of the band structure. This has also 
been suggested in Ref. [0, |2(| to explain the large inverse isotope effect observed in 
Uranium. In this context, it would be interesting to determine the uranium isotope 
effect in heavy-fermion systems (one should note, however, that there is only one report 
of the isotope effect in uranium flT2f). 

Another effect of the band structure arises if one assumes that the electronic density 
of states varies strongly on a scale given by Q, the BCS energy cutoff (e.g., in the 
presence of a van Hove singularity). To obtain the BCS expression (^) we assumed 
that the electronic density of states (DOS) is constant in the energy intervall [— Q, Q] 
around the Fermi energy. A is then given by N(Ep)V, where N is the DOS at the 
Fermi level and V is the attractive part of the electron-electron effective interaction. In 
a more general case, however, one has to consider the energy dependence of the DOS. 
Though a purely electronic parameter, the energy dependence of the DOS appears to 
influence the isotope effect [f2~T|1 . One can understand this effect qualitatively within a 



crude extension of the BCS model. Instead of taking the DOS at the Fermi level, one 
replaces N(Ep) by < N(e) >n, the average value of the DOS N(e) over the interval 
[— O, fl] around Ep. Obviously, this average depends on the cutoff energy Q of the 
pairing interaction. Isotopic substitution modifies Q which affects < N(e) >q, and 
thus \=< N > n V and T c . 

A more general analysis of this effect within Eliashberg's theory was given in 



Ref. ||21|| . They find that whereas the IC can reach values above 0.5, its minimal 
value obtained for reasonable choices of the parameters never reaches the small values 
(~ 0.02) found in optimally doped high-T c compounds. However, more recent stud- 
ies of the influence of van Hove singularities on the isotope coefficient show that such 
small values can be obtained, because in this scenario the cutoff energy for the effective 
interaction between charge- carriers is given by min(Ep — E v h;Q) (where E v h is the 
energy of the van Hove singularity and Q is a characteristic phonon energy) f22[ , [23| , p4|j . 
If the cutoff energy is given by Ep — E v h and is thus electronic in origin, T c displays 
no shift upon isotopic substitution. Taking into account the Coulomb repulsion it has 
been shown[24] that one can even obtain a negative IC. 

Another way to extend the results presented earlier is to consider an anisotropic 
Eliashberg coupling function a 2 F(q, uj). Generally it is assumed that the system is 
isotropic and the coupling function can be averaged over the Fermi surface, leading to 
a q-independent Eliashberg function. If the system is strongly anisotropic (as is the case 
of high-T c superconductors), the average may not lead to an accurate description of the 
situation. The isotope effect has been studied for different q-dependent form factors 



entering a F in Ref. f25 |. They show that a small isotope coefficient as observed in 



YBa2Cu307_«5 can be obtained for anisotropic systems. Nevertheless, to obtain high 
critical temperatures at the same time, they are forced to introduced another, electronic 
pairing mechanism. Such a situation is further discussed below (see Sec. |2.6|) . 

Note that the theories presented in this section give only a qualitative picture of 
the situation encountered in high-temperature superconductors and do presently not 
account for the isotope effect observed in various systems (see also Refs. JTL], f26[|). 



2.4 Polyatomic Systems 



Until now we focused mainly on the study of monoatomic systems in which the attrac- 
tive interaction leading to the formation of pairs is mediated by phonons. All previous 
effects are naturally also encountered in polyatomic systems. However, the presence of 
two or more elements in the composition of a superconductor has several direct conse- 
quences for the isotope effect that we summarize in the following. The first consequence 
of a polyatomic system is that the characteristic frequency Q which determines the value 
of T c depends on the mass of the different ions involved Q = Q,{M\, M 2 , ■ ■ •). Obviously, 
the dependence Q ~ M~ ar for element r = 1, 2, ... must not be equal to a r = 0.5. It is 
thus not surprising if the partial isotope coefficient (obtained by substituting one type 
of ions for its isotope; see Eq. (§)) differs from the textbook value 0.5. 

Let us consider the example of a cubic lattice with alternating masses M\ and 



M 2 []27|, p8fl . The acoustic and optical branches can be calculated analytically [^7]]. From 



this one obtains the following partial isotope coefficients: 

ai = ^ { K± R J" +2R i M >} , (12) 
1 (f!=F/?M+) 

°> = 2mr,\ K± to T /?M 5 ) > (13) 




where = Mf 1 - M 2 _1 , Mjtj = Mf 1 + M 2 _1 . The upper sign stands for the Debye 
frequency and the lower sign has to be considered when the characteristic frequency 
is given by the optical phonon. K = J2 x , y ,z Ki is the sum °f the force constants and 
Kl = J2x, y ,z Ki cos(gL) (q — for the optical branch and q = n/L for acoustic mode; 
L is the lattice constant). One notes that the total isotope coefficient is given by 
c^tot = cut + ci2 = 0.5 whether one takes the acoustical or the optical phonon as the 
characteristic energy for the determination of T c . 

Given the previous result, one can ask if there exists a maximal value of the IC. It is 
often stated in the literature that 0.5 is the maximal value that the isotope coefficient 
can reach within a harmonic phonon, and electron-phonon induced pairing model. Al- 
though to the knowledge of the authors none of the systems studied so far seems to 
contradict this assertion, one should note that there is no proof of this statement. It 



was shown in Ref. |29| that for a polyatomic system atot = H T o> r = 0.5 (defined in 
Eq. (p|)) if one assumes that fi* = and all masses of the unit cell are subject to the 
same isotopic shift (i.e., AM r /M r takes the same constant value for each element r of 
the system, see Eq. (0)). These assumptions hold approximately for certain systems 



as, e.g., the Chevrel-phase MoeSes material |29|1, but are certainly not valid for example 



in high-temperature superconductors. It is thus not clear if the total isotope coeffi- 
cient, Eq. (0), is indeed always 0.5. The general proof of such a statement requires the 
knowledge of the polarization vectors and their derivatives with respect to the isotopic 
masses, both of which have to be calculated for each specific system studied. 

Another important remark concerns the value of the (very small) IC for optimally 
doped high-T c superconductors. Experimentally one measures generally the partial 
oxygen isotope effect (in some cases also the Cu and Ba IE). Since the unit cell of a 
high-temperature superconductor contains many different atoms, one expects values of 
the partial IC that are significantly smaller than 0.5. A crude estimate can be given 
by observing that under the assumption of a same contribution of each atom to the 
isotope coefficient (which is certainly not the mentioned above) the value of the 

isotope effect for oxygen in YBCO would be a ~ 0.5/iV w 0.04, where iV is the number 



of atoms per unit cells (N = 13 in the case of YBCO). Since there are 6 — 7 oxygen 
atoms in the unit cell, one should multiply a ox by this number. 

If, in addition to the multiatomic structure of the system, one takes into account 
the effect of Coulomb interactions and/or the fact that the charge- carriers may strongly 
couple to certain phonons only one can obtain values of the order observed in high- 
temperature oxides. The second possibility has for example been studied in the context 
of the Chevrel-phase compound Mo 6 Se 8 [^| or organic superconductors |JD|. It has also 



been studied for high-T c superconductors in Ref. |3T], j32j1 - They conclude that the 
oxygen isotope effect can be well described for reasonable parameters in the case of 
La 2 _ 3; Sr 3 ;Cu04, but that unphysical values of the coupling have to be considered to 
describe the partial oxygen isotope effect of YBa 2 Cu 3 07_5. A careful study of the 
effect of Coulomb interactions on the IE of high-T c superconductors has not been done 
yet. 

On thus concludes from the previous considerations that it is not impossible to 
explain the small values of the IC observed in conjuction with high T c 's in high- 
temperature materials, if one considers the combined effect of a polyatomic system, 
the Coulomb interaction and the band structure. 



2.5 Anharmonicity 



Anharmonic effects play an important role in materials such as PdH(D)j33], [34 



Mo6Se 8 [^9| or, according to some theories p6|, p5| , [36], [37|], in high-temperature supercon- 
ductors. Among other consequences the presence of anharmonicities affects the value of 
the isotope coefficient. We present here two different aspects of the anharmonic isotope 
effect: anharmonicity of the characteristic phonon mode (for PdH) and volume effects 
(for Molybdenum). 

Quantities that are independent of ionic masses in the harmonic approximation, 
can become mass dependent in the presence of anharmonicity. We mention three such 
properties that are of interest for the isotope effect in superconductors: the lattice force 
the unit-cell volume ||3l 



constants fC;[|34 
A 



39 1 and the electron-phonon coupling function 
6| |35[. In the first case, it was shown that the account of Coulomb interactions alone 
see above) cannot explain the inverse (a < 0) hydrogen isotope effect observed in PdH 
see Table 1). On the other hand, a change of 20% of the lattice force constants when 



replacing H for Deuterium (D) in PdH was infered from neutron-scattering data ||34|| . 
This results from large zero-point motion of H as compared to D. Taking this fact into 
account in the Eliashberg formalism allows one to obtain quantitative agreement with 



the experiment [|34 ] . 

Another effect due to the presence of anharmonicity is the isotopic volume 



effect |8], |39 



This effect was discussed in the context of Molybdene and is related 
to the difference in zero-point motion of the two isotopes. In such a system one can 
write the IC as a = a bcs + a voi with 



M AV AT C 



(14) 



where B = —VdP/dV is the bulk modulus, P is the pressure and AV — V* — V is 
the volume difference induced by isotopic substitution. In the case of Molybdenum 
the volume isotope effect amounts to a vo i ~ 0.09 and accounts for ~ 27% of the total 
isotope coefficient, which is not a negligible effect. 

The presence of a volume isotope effect has also been suggested in the case of PdH, 
Ceo materials as well as in Pb|p8||. 



The other major effect of anharmonicity is the appearence of an ionic-mass de- 
pendent electro n-phonon coupling function A|26|, |35], [42| . In this case not only the 
prefactor Q of Eq. ([|), but also the exponent depends on ionic masses. As stressed 
earlier in the context of band structure effects, this fact can lead to a strong deviation 
from olbcs — 0.5, even for monoatomic systems. An explicit expression for the motion 



of the oxygen in a simple double-well potential was obtained in Ref. [3l| for high-T ( 



superconductors. The model was extended further in Refs. |35[b and [26 and it was 
shown that the isotope coefficient can be much smaller than 0.5 with moderate values 
of the coupling constant or exceed 0.5 for (very) strong couplings (the function a(A) 
goes through a minimum at A ~ 1; see Ref. f35|). Regarding high-T c superconductors, 
the anharmonic model could explain qualitatively the behaviour of the isotope coef- 
ficient if one considers the fact that the coupling function A also depends on doping 
(see Ref. ||3| for this dependency). Assuming that these systems are strong-coupling 
superconductors, and that the coupling decreases upon doping, one can describe qual- 
itatively the behaviour of a(T c ) for underdoped systems. However, it is difficult to 
obtain quantitative agreement especially in the optimally doped and strongly under- 
doped regime. In the first case the coupling has to be intermediate to obtain small a's 
but then T c is also small. In the second, strongly underdoped case, the IC can exceed 
0.5, but the coupling has to be very strong. One way to solve this problem is to assume 
that the superconducting pairing is mediated by an additional, non-phononic channel 
(see Ref. ||2[ and below). 



The influence of anharmonicity on the electron-phonon coupling A and on the hop- 
ping parameter (for a tight-binding type of lattice) has also been considered in the 
context of 12 C ^ 13 C isotope substitution in A 3 Cqo materials (A =Na, Ru; see, e.g., 
Ref. |^0| and references therein). Experimentally the isotope coefficient was shown to 



vary between 0.37 to 1.4, depending on the isotopic substitution process (in the first 
case each Cqo molecule contains an equal amount of 12 C and 13 C isotopes whereas in 
the second case the system contains Ceo-molecules that are composed of either pure 
12 C or pure 13 C atoms) @. 



2.6 Non-Phonon and Mixed Mechanisms 

The superconducting transition can also be caused by a non-phonon mechanism. His- 
torically, the introduction of the electronic mechanisms started the race for higher 



T c 's. The pairing can be provided by the exchange of excitations such as excitons, 
plasmons, magnons etc... (see, e.g., Ref. []25]]). m general, one can have a combined 
mechanism. An electronic channel can provide for an additional contribution to the 
pairing, as is the case, e.g., in a phonon-plasmon mechanism |44[ . In the context of the 



isotope effect, the mixed mechanism may provide for an explanation for the unusual 
occurence of high T c 's and small isotope coefficients a (see above). 

Let us initially consider the case where only electronic excitations mediate the pair- 
ing interaction and the Eliashberg function a 2 F can be approximated by a single peak 
at energy Q e , the characteristic electronic energy. The theory considered in weak cou- 
pling yields then a relation of the form T c ~ f2 e e -1/,A . For an electronic mechanism VL e 
is independent of the ionic mass. Therefore, the isotope coefficient of T c vanishes for 
such cases. 

If, on the other hand, the superconducting state is due to the combination of phonon 



and a high-energy excitations (as, e.g., plasmons [|3]) then one has to include both 
excitations in a 2 F. In this case, the simplest model considers an Eliashberg function 
composed of two peaks, one at low energies (Qo) for the phonons and one at high 
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energies (f^) for the electronic mechanism. For such a model, one obtains the following 
expression for T c ||45|: 



T c = l.UQM 1 exp (_M^T| (is) 
\ Po + pi J 

where h is a slow varying function of pi and fli, fi = pi/ (po+pi) and pi = Ai/ (1+Ao+Ai) 
{i = 1, 2). In weak coupling, h ~ 1 and the exponent is independent of the ionic masses. 
The resulting isotope coefficient for this mixed mechanism takes then the form: 

If only the phonon mechanism is active Ai = and one recovers the BCS value 0.5. 
Otherwise, the presence of the non-phononic mechanism reduces the value of the isotope 
coefficient, while enhancing T c . For Ao <C Ai one has a ~ 0. A joint mechanism 
would thus allow one to explain the small values of the IC for optimally doped high-T c 
superconductors since the coupling to non-phononic degrees of freedom would provide 
for high T c 's while it would reduce, together with the Coulomb repulsive term, the value 
of the IC. 

Another way to include non-phononic contributions to the pairing mechanism, is 
to consider a negative effective Coulomb term p* in the Eliashberg theory |28|, [46| , pi]] . 
The pseudo-potential p* can thus be seen as an effective attractive pairing potential 
contribution and thus supports superconductivity. Such a negative p* is for example 
obtained in the plasmon model |44|, |28|, |46| . One notes that the small values of the 
partial isotope coefficient in optimally doped high-T c materials can easily be explained 
within such a model. 

Several other non-phononic as well as combined mechanisms have been proposed 



[pT] , |42| , [17], f|8| that also lead to a reduction of the isotope effect. We refer to the 



literature for the details. 



2.7 Isotope Effect of Properties other than T c 

All previous considerations were concerned with the isotope effect of the superconduct- 
ing critical temperature T c . One can ask if there are other properties displaying an 
isotope effect in conventional superconductors. Naturally, every quantity depending 
directly on phonon frequencies will display such an effect. Let us consider here a prop- 
erty that will be studied further in the next sections, namely the penetration depth of 
a magnetic field 5. In the weak-coupling London limit the penetration depth is given 
by the well-known relation 

2 = mc 2 = mc 2 

4nn s e 2 4nn<p(T/T c )e 2 ' 1 ' 

where m is the effective mass. n s is the superconducting density of charge carriers, 
related to the normal density n through n s = n<p(T/T c ). The function (p(T/T c ) is a 
universal function of T/T c . For example, (p ~ 1 — (T/T c ) 4 near T c , whereas (p ~ 1 near 
T = (in the absence of magnetic impurities; their influence is discussed in section |3|). 

Eq. ( |TTD does not depend explicitely on phonon frequencies. Nevertheless, it can 
display an isotopic dependence through T c . In fact, this dependency is common to all 
superconducting properties within the BCS theory. Indeed, all major quantities such 
as heat capacity, penetration depth, critical field, thermal conductivity etc... can be 



expressed as universal functions of T c . As a result, all these quantities display a trivial 
isotopic shift, caused by the isotopic shift in T c . The value of the shift is growing as one 
approaches T c . Furthermore, the shift vanishes for T —>■ 0. As will be shown in section 
[| and following the situation is very different in high-temperature superconductors and 
manganites. For example, the influence of non-adiabaticity on charge-transfer processes 
results in an unconventional ionic mass dependence of the charge-carrier concentration 
n = n(M), which, according to Eq. (|TTD, leads to a new isotope effect (see sec. ||). 

So far, we have presented various theories explaining the deviation of the measured 
isotope coefficient in conventional superconductors from the value a = 0.5 derived 
within the BCS model. The remainder of this review is devoted to the study of three 
factors that are not related to the pairing interaction but affect the isotope coefficient: 
magnetic impurities, proximity contacts, and non-adiabatic charge-transfer processes. 
We show that they can strongly modify the value of the isotope coefficient and can even 
induce an isotopic shift of superconducting properties such as the penetration depth 
5. In the last part of this review, we then apply the theory to the case of the oxygen 
isotopic substitution in high-temperature superconductors. 



3 MAGNETIC IMPURITIES AND 
THE ISOTOPE EFFECT 

The presence of magnetic impurities strongly affects various properties of a supercon- 
ductor. It has been shown [pED|, pTi| ] that because of magnetic impurity spin- flip scattering 
processes, Cooper pairs are broken and thus removed from the superconducting conden- 
sate. Several properties are affected by the pair-breaking effect. Magnetic scattering 
leads to a decrease of the critical temperature T c , the energy gap, the jump in the 
specific heat at T c , etc... At some critical impurity concentration tim = nu,cr super- 
conductivity is totally suppressed. Moreover, there exists an impurity concentration 



71m = n-M,g < nM, cr (for conventional superconductors um, 9 = 0.9nM,cr, see Ref. [|9|, |50"H ) 
beyond which the superconducting state is gapless. Pair-breaking also leads to an in- 
crease of the penetration depth S. We show in the following sections that the presence 
of magnetic impurities leads to a change of the isotope coefficient of T c and S. 



3.1 The Critical Temperature 

The change of T c induced by magnetic impurities is described in the weak-coupling 
regime by the Abrikosov-Gor'kov equation J49|j: 

1*^=^ + 7^-^ , (18) 



TJ r V2 '7 V2 

where 7 S = T s /2irT c and T c (T c0 ) is the superconducting critical temperature in the 
presence (absence) of magnetic impurities. Y s = T s n M is the spin-flip scattering ampli- 
tude and is proportional to the magnetic impurity concentration i%m (r s is a constant). 

It is easy to derive a relation between the isotope coefficient «o in the absence of 
magnetic impurities [a = — (M/AM)(AT c0 /T c0 ), Eq. (|TJ)] and its value in the pres- 
ence of magnetic impurities [a = —(M/AM)(AT C /T C )]. From Eqs. (|l|) and (|18|) one 
obtains |5T], f§ § 0: 



am W(% + 1/2)7 S ' (19) 
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where ip' is the derivative of the psi function and is a positive monotonous decreasing 
function of j s . The fact that ao and a are not identical relies on the essential feature 
that the relation between T c q and T c (Eq. ([18])) is non-linear. It is interesting to note 
that magnetic scattering decreases T c (reduction of the condensate) but increases the 
value of the IC (a m > «o; see Eq. (|19D). 

An important consequence of Eqs. fll8f) and (|i~9[) is that one can write the IC as a 
universal function a(T c ). This relation does not contain any adjustable parameter, but 
depends solely on the measurable quantities a and T c . The universal curve is shown 
in Fig. 2. Note also that a m can, in principle, exceed the value ascs — 0.5. Thus, the 
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Figure 2. Isotope coefficient of T c in the presence of magnetic impurities. Solid line: 
Universal dependence a(T c ) (normalized to ao; see Ref. ]|]); Points: Isotope effect for 
YBa2(Cui_ a; Zn :r )307 obtained with various experimental techniques and normalized to 
a = 0.025 (from Ref. ||, g). 

observation of large values of a are not necessarily related to the pairing mechanism. 



3.2 Zn-doped YBa 2 Cu 3 
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Zn substitution for Cu in the Cu02 planes of YBCO has attracted a lot of interest, since 
it leads to a drastic decrease in T c and a strong increase of the penetration depth 
|54j (superconductivity is destroyed with ~ 10% of Zn). In addition, this decrease of T c 
is accompanied by an increase of the isotope coefficient (see Ref. [53, Q and Fig. 2). We 
think that the peculiar behaviour of Zn doping is related to the pair-breaking effect. 
It has been established by several methods (see Refs. [^3|, [54]]) that Zn substitution 
leads to the formation of local magnetic moments (~ 0.63/i^/Zn, where fis is the Bohr 
magneton) in the vicinity of Zn. 

Fig. 2 displays the experimental results obtained for the oxygen isotope coefficient 
of Zn-doped YBCO[^, U (normalized to cto — 0.025|57]]). One can see from this figure 
that the agreement between the theoretical dependence a(T c ) and the experiment is 
very good. Note that the uncertainty in the data is growing as T c — > (see Ref. |)6| ) . 
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3.3 The Penetration Depth 



The isotope effect of the penetration depth is more complicated than the isotope effect 
of T c , since it appears to be temperature dependent. Here we present the main results 
near T c and at T = and refer the reader to Ref. [§, |6| for details. Near T c , the 
penetration depth is given in second order of the superconducting order parameter 
A = A(T,r s )byH: 

5- 2 = a^-a^ls + \) , (20) 

where a = 4:(Jn/c (oat is the normal state conductivity), ((z, q) = J2 n >o V( n + c l) z an d 
7 S was defined in Eq. (0). One can calculate the isotope effect from Eqs. (|3|), (p0|) and 
the analytical expression for A near T c (see Refs. |5(| 



a z = 2nd -r 



2 ... N l-C 2 + (l-r) I-C2 + C 



,2^1 



C3-C4 S A 



= 2rj-i (21) 



with C 2 = 7r 1 C(^ ) 7, + 1/2), 2 = 1, 2, ... and r = T/T c . 

As mentioned in Sec. 0, 5 experiences a trivial BCS isotope effect near T c . Since we 
are only interested in the unconventional IE resulting from the presence of magnetic 
impurities, we substract the BCS isotope effect by calculating the isotope coefficient 
P m of 5(T, T S )/5(T, 0). Near T c this coefficient can be written as /3 m = (3 m — /3q, where 
Pm (Po) is the isotope coefficient of S in the presence (absence) of magnetic impurities. 
One can showH that (near T c ) f3 m can be rewritten in the form: 



where 



P m = (R x - R )a m , (22) 



Ri = -I/ 1 = -lf^-^ + 2$L-l) , (23) 

= - \— h = ~[l ~ ^(7. + 1/2)7.] /o , (24) 
2 a m 2 



and /o = (3 — r 2 )/(l — r)(3 — r). The functions N\, D% are defined in Eq. (|2l| ) and 

iV 2 = 3(1 - r) 2 [C2 - 2C 3 + C4I + r(2 - r) - C2 
^2 = 2 (30 - Cs - 2Cs 



This relation is valid near T c (where A is small) and for impurity concentrations such 
that AT C /T C < 1. 

One notes first that the isotope coefficient of 6 is proportional to the IC of T c . As 
for the isotope coefficient of T c (Eq. (fl9f) ) all quantities can either be obtained from 
experiment (e.g., a , T c0 and T c ) or calculated self-consistently using Eqs. ([Df) and 
(|l9|). There is thus no free parameter in the determination of (3 m . Fig. 3 displays the 
IC Pm (normalized to ato) as a function of T c (i.e. on the concentration um) for fixed 
values of the temperature. Setting a = 0.025, one obtains the relation a(T c ) expected 
for YBCZnO. It would be interesting to measure the isotope effect of S near T c in this 
system, to verify our predictions. 

An interesting property of Eq. fl2"2] ) that is shown in Fig. 3 is the fact that the 
isotope coefficient of 5 is temperature dependent; an unusual feature in the context 
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Figure 3. Isotope coefficient j3 m near T c (normalized to ato) as a function of T c (i.e. of mag- 
netic impurity concentration) for T/T c = 0.75 (solid line), 0.85 (dashed) and 0.95 (dotted). 
Multiplying by oq = 0.025 gives the IC expected for YE^^ui-^Zn^C^^. 

of the isotope effect. We will see in Sec. f| that a temperature dependent IC is also 
observed for a proximity system. 

As can be seen by comparing Figs. 2 and 3, the qualitative behaviour of /3 m (T c ) is 
similar to a m (T c ) [Eq. (|i~9D with ao = a p h] but with opposite sign. Furthermore, in the 
absence of magnetic impurities one has $ m — 0, whereas a m (T c ) = a p h- 

Note finally, that since R\ — Ro < 0, the IC of the penetration depth is always 
negatif when magnetic impurities are added to the system. This conclusion might 
not hold in certain cases if the non-adiabatic channel is also included (then one has 
Ptot = Pm ~ Po + Pna and (3 na is positive, see Sec. g). 

All previous considerations have been done near T c . At T = 0, the BCS contribution 
arising from the T c dependency of the superconducting charge-carrier density n s (see 
Eq. ( PD ) vanishes, because (f(T = 0) = 1. The only contribution to the IE is thus 
due to magnetic impurities. In the framework of the Abrikosov-Gor'kov theory one can 
write the isotope coefficient at T = in the form: 

P m = R a , (25) 

where Ro is given in appendix. Note that Ro is a negative function of T 2 (the direct 
scattering amplitude) and T s = Fx — T 2 (the spin-flip, or exchange scattering amplitude) 
as defined by Abrikosov and Gor'kov||49|. The IC of 5 is thus negative both near T c and 
at T = 0. This has to be seen in contrast to the IC of T c which is always positive (see 
Fig. 2). Fig. 4 shows the universal relation P(T C ) for two values of F 2 . One remembers 
that T c is determined by the magnetic impurity concentration. 

There are two major differences between the results obtained near T c (Fig. 3) and 
at T = (Fig. 4). First, the direct scattering amplitude F 2 appears only at T = 0. 
Secondly, the IC near T c is proportional to a m given by Eq. fllDQ, whereas the IC at 
T = is function of ao, the IC of T c0 in the absence of magnetic impurities. This 
difference is due to the fact that ip = 1 at T = (see Eq. (|T7|) ) and the penetration 
depth does consequently not depend on T c (the critical temperature in the presence of 
magnetic impurities). 
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Figure 4. Isotope coefficient of the penetration depth at T = in the presence of magnetic 
impurities (normalized to «o)- There is no free parameter in the relation f3 m (T c ). Solid and 
dashed lines are for Tz/Tg = 10 and 50 respectively (see Ref. @). 



As shown in Fig. 2, the IC of T c fits the data of Zn-doped YBCO. On the other hand, 
there are no experimental data for conventional superconductors. Furthermore, the 
magnetic impurity contribution to the IE of 8 shown in Fig. 3 has never been measured. 
It would thus be interesting to perform measurements of these effects, especially since 
they can be described by universal relations. 



4 ISOTOPE EFFECT IN A 
PROXIMITY SYSTEM 

In this section we consider another case, a proximity system, in which a factor not 
related to lattice dynamics induces a change of the isotope coefficient of T c and 8. 
Consider an S — N sandwich where S (N) is a superconducting (normal) film. The 
value of the critical temperature T c and the penetration depth 8 differ substantially 
from the values T c0 and 8$ of the isolated superconductor S^ty, |59|j . We show in the 
following that the proximity effect also influences the isotope coefficient of T c and S. 



4.1 The Critical Temperature 



Let us consider a proximity system composed of a weak-coupling superconductor S of 
thickness L$ and a metal or a semiconductor N of thickness (e.g. Nb-Ag). In the 
framework of the McMillan tunneling model[]T7|], which can be used when 8 < 
(6v = hv F . N /2nT is the coherence length of the iV-film as defined in Ref. [|SC(] ), the 
critical temperature T c of the S — N system is related to the critical temperature T c o 
of 5" by the relation [Q: 



vrT, 



c0 



■ cO 



2^u 



P 



v N L 



N 



(26) 



where 7 ~ 0.577 is Euler's constant. The value of u is determined by the interplay of 
the McMillan tunneling parameter r = r,s7v + rArs and the average phonon frequency 



Q. For an almost ideal S — N contact (T 3> Q) one has u ~ Q. In the opposite limit, 
when r -C Q one obtains u ~ T (see Ref . |58| ) . 

Let us first consider the case r ^> Q. In the BCS model, one has T c0 oc fl Thus, 
according to Eq. (^) T c and T c o have the same dependency on the ionic mass in this 
limit (because T c q/u is then independent of Q). This results in the simple relation 
otprox = «o for the IC of T c («o is the isotope coefficient of T c for the isolated S'-film). 

The opposite situation where r <C Q is more interesting since u = T is independent 
of ionic masses. Using Eqs. ([!]) and ([26|) one obtains for the isotope coefficient of T c in 
the limit r < Q: 

Otprox = «0 ^1 + ■ ( 27 ) 

One notes first that whereas the presence of a normal film on the superconductor 
decreases T c , the film induces an increase of the isotope coefficient of T c . The same 
is true for the presence of magnetic impurities and for the non-adiabatic IE when 
dT c /dn > (see next section). The second interesting feature of Eq. ( p7|) is that one 
can modify the value of a prox by changing the thicknesses of the films. For example, if 
vn/^s = 0.8 and L^/Ls = 0.5 then a prox = 0.28. By increasing the thickness of the 
normal film such that Ln = Ls, one obtains a prox = 0.36. 

We stress the fact that, as in the previous section on magnetic impurities, the change 
of the IC is due to a factor not related to lattice dynamics. Therefore, there is no reason 
for a prox to be limited to values below 0.5. 

To the best of our knowledge the change in the IC caused by the proximity effect has 
never been measured, even in conventional superconductors. It would be interesting to 
carry out such experiments in order to observe this phenomenon. 



4.2 The Penetration Depth 

The proximity effect also affects the shielding of a magnetic field. The most dramatic 
effect of the normal layer on the penetration depth is seen in the low-temperature regime 
(T/T c < 0.3). Although the penetration depth of a pure conventional superconductor 
is only weakly temperature dependent in this regime (5~ 2 ~ if = 1 — T 4 /T c 4 ) the 
presence of the normal layer induces a temperature dependence through the proximity 
effect. For the same S — N proximity system as considered above the penetration depth 
is given by [|5"9|l 

5- 3 = a N $ , (28) 



where a at is a constant depending only on the material properties of the normal film 
(it is ionic-mass independent) and 



$ = ^Ezi-T-7> Pn = l + et^xl + l . (29) 

n>0 ^nPn ~r -L 

x n = u n /es(T) with u n = {2n + \)ixT (the Mastubara frequencies) and e s {T) is the 
superconducting energy gap of S. In the weak-coupling limit considered here, es(0) = 
£7tT c0 with e ~ 0.56. The dimensionless parameter t = £/So with £ = Ln/L q and 
So = T/7tT c o (r ~ 1/Lq is the McMillan parameter [|T7|1 ). Ln and L$ are the thickness 



of the normal film and some arbitrary thickness, respectively (in the following we take 
L = L s , the thickness of the superconducting film). From Eqs. (131) and (E3) one 



I 1-7 



obtains the IC of the penetration depth for the proximity system: 



2a 



prox 



3$ 



n>0 



/v>2 ~,2 



et 



XnPi + 1 



PnV^n + l , 



(30) 



where «o is again the IC of T c0 for the superconducting film S alone. This result has 
three interesting features. First, as for the case of magnetic impurities, the IC of T c 
and 5 have opposite signs. This observation is valid as long as one does not mix the 
different channels presented in this work (magnetic impurities, proximity effect and 
non-adiabaticity). The addition of a non-adiabatic contribution (next section) may 
lead to a different conclusion. Secondly, one can see from Eq. (|30"D that (3 prox depends 
on the proximity parameter t. One can thus modify the IC either by changing the ratio 
I = Ln/Ls or the McMillan tunneling parameter T (e.g. by changing the quality of the 
interface; see Ref. 0). Fig. 5 shows this dependence for different temperatures. One 
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Figure 5. Isotope coefficient (normalized to ao) of the penetration depth for a proximity 
system as a function of t, for T/T c = 0.1 (solid), 0.2 (dashed), 0.3 (dotted). See Ref. ||. 

notes that, contrary to the value of a prox , the IC of the penetration depth decreases with 
increasing ratio t. Finally, one notes that the IC (3 prox is temperature-dependent. This 
unconventional feature was also observed for the IC of 5 in the presence of magnetic 
impurities near T c . Fig. 6 shows the temperature dependence for different values of the 
parameters. The trend is similar to the case of magnetic impurities: |/3 pr . ox | increases 
with increasing T. Note, however, that the two effects are calculated in very different 
temperature ranges (near T c and near T = 0). A complete description of this effect can 
be found in Ref. H. 



5 NON-ADIABATIC ISOTOPE EFFECT 



The non-adiabatic isotope effect introduced in Ref. was used in Refs. [0, f|, ||, [F| 
to describe the unusual behavior of the isotope coefficient in several high-temperature 
superconductors. The theory also allowed us to describe the large isotopic shift of the 
ferromagnetic phase transition in manganites [0 . The concept of the non-adiabatic IE 
relies on the fact that in the above-mentioned systems charge-transfer processes involve 
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Figure 6. Isotope coeff. (3 prox (normalized to olq) for a proximity system as a function of 
T/T c for S" = 0.2: I = 1 (solid) I = 0.5 (dashed) and S = 1: 1 = 1 (dash-dotted) I = 0.5 
(dotted). See Ref. @. 



a non-adiabatic channel. In other words, the electronic ground state of the group of 
atoms over which the charge-transfer takes place is degenerate, leading to a dynamic 
Jahn- Teller (JT) effect (E JT < hu, where E JT is the JT energy; see, e.g., Ref. f61|j ). 



The concept of the non-adiabatic isotope effect can be understood by considering 
two examples: the apex oxygen in high-T c oxides and the motion of charge- carriers in 
manganites. Let us first consider high-T c materials. This system can be seen as a stack 
of charge reservoirs (CuO chains) and conducting subsystems (Cu0 2 planes). Charges 
are transfered from the reservoir to the conducting layers through the apex oxygen 
that bridges the two subsystems. Several experiments have shown [|6^ , |63| , |64| that 



the apex oxygen displays a non-adiabatic behaviour. The ion oscillates between two 
close positions that correspond to two configurational minima of the potential energy 
surface. These two minima are due to the Jahn- Teller effect since we are dealing with 
the crossing of electronic terms (see Fig. 7). 

The Jahn- Teller effect leads to a "double-well" type potential (which should not 
be confused with the double-well appearing when the crystal is anharmonic; here we 
consider energy terms crossing and for each term we use the harmonic approximation). 
This structure has a strong impact on charge-transfer processes in these materials. 
Indeed, the motion of charges from the chains to the planes occurs through the apex 
oxygen. Thus, the charge-transfer process involves the motion of the non-adiabatic ion. 
The immediate consequence of this observation is that the density of charge- carriers n in 
the Cu0 2 -planes (the conducting subsystem) depends on the mass of the non-adiabatic 
ions involved n = n(M) (see below). 

The situation encountered in manganites is different from the example just de- 
scribed. In these materials, there is no transfer between a reservoir and a conducting 
subsytem. Instead, the motion of charge-carriers occurs on the Mn-O-Mn complex 
that displays a non-adiabatic behaviour |65[ . As a consequence, the electron hopping 



involves the motion of the ions and depends on their mass. It was shown in Ref. 
that this fact can account for the unexpectedly low ferromagnetic critical temperature 
T c j |66|1 (which is determined by the hopping of the electrons between Mn ions) as well 
as for the huge isotopic shift of T c /|6l 




Figure 7. Potential energy surface (PES) as a function of the configurational coordinate R. 
The crossing of electronic terms leads to a dynamical Jahn- Teller distortion (two minima of 
the PES). 



Note that although the charge-transfer processes are different in high-temperature 
superconductors and in manganites, the formalism described below and in Ref. [p8 



applies equally well to the two systems. In the following we refer to the example of 
the apex oxygen. The case of manganites can be mapped by replacing "reservoir" and 
"conducting subsystem" by "Mn ions" . 

It is important to realize that the dependence n(M) is specific to systems where 
the motion of charges occurs through non-adiabatic ions. In metals, the charge- carriers 
are not affected by the motion of the ions because of the validity of the adiabatic 
approximation. We also stress the fact that the apex oxygen is given as a simple 
example but, as will be shown below, the theory is not limited to this case. Finally, one 
should note that a phenomenological theory of the IC based on the assumption that 
n = n(M) has been proposed in Ref. . 



To demonstrate that in the presence of non-adiabatic charge-transfer one has 
n = n(M), one best uses the so-called diabatic representation [|70|1 (see also Ref. p8|). 



Here we only recall the main steps of the calculation. Let us assume that, because of the 
degeneracy of electronic states (or Jahn- Teller crossing), the potential energy surface 
E(R) (R is the relevant configurational coordinate) of the group of non-adiabatic ions 
(such as the apex O, in-plane Cu and O for high-T c materials or Mn and O ions in man- 
ganites) has two close minima (see Fig. 7). In the diabatic representation the total wave- 
function \l/(r, R, t) (r is the electronic coordinate) is written as a linear combination of 
the wave-functions in the two crossing potential surfaces. The total wave-function can- 
not be decomposed as a product of electronic and lattice wave-functions. On the other 
hand, one can write the total wave-function as a sum of symmetric and antisymmetric 
terms. The energy splitting between these symmetric and antisymmetric terms, which 
corresponds to the inverse lifetime of oscillations between the configurations (minima 
of the potential energy surface), has the form|| ||, [7|, |68fl : 

i/ 12 =<^ 1 |i/ e |^ 2 >~L F 12 , (31) 



on 



where H e is the electronic part of the total Hamiltonian and 



L = L(Ro) = |dr^*(r,R)# e ^ 2 (r,R) , (32) 
F12 = J rfR$t(R)$ 2 (R) • (33) 

The last equality in Eq. fl3~T[ ) is obtained under the assumption that the electronic 
wave-function ipi (i — 1, 2) is a slowly varying functions of R. L Q can then be evaluated 
at R , the crossing of electronic terms, and taken out of the integral over R. It does 
not depend on ionic masses. On the other hand, the important Franck-Condon factor 
(|33| ) depends on the lattice wave- functions $i(R), and thus on ionic masses. Given the 
perturbation, Eq. ([31]), one can calculate the probability of finding the charge- carrier 
in the conducting layer || |, [7], |68| . Note that non adiabaticity affects the bandwidth 
of the reservoir and conducting bands rather than the chemical potential (which is the 
same because of thermodynamic equilibrium). Thus, in the case of YBCO, the plane 
and chain bands have the same chemical potential, but different widths. This affects 
the charge- carrier concentration in the two subsystems. 

Qualitatively, the charge transfer can be visualized as a multi-step process. First, the 
charge carrier can move from the reservoir to the group of ions (e.g. from the chains to 
the apical oxygen in YBCO). Then, the non-adiabatic ions tunnel to the other electronic 
term (\l/i — > ^2)- As a final step, the charge carrier can hop to the conducting layer 
(from the apical oxygen to the Cu02 planes). The crucial point of the theory is that 
as a result of this multi-step process, the probability P of finding the charge-carrier in 
the conducting subsystem depends on Eq. ( pi]) and thus on the Franck-Condon factor 
F12 = -Fi 2 (M), Eq. (|3~3"1), that depends on ionic masses. In other words the charge- 
carrier concentration n, which is proportional to the probability P, depends on the 
ionic mass M. It is this unusual dependence n = n(M), found in high-temperature 
superconductors and manganites, that is responsible for the unconventional isotope 
effect of T c and S. 



5.1 The Critical Temperature 

Given that n = n(M), we can write the isotope coefficient of T c as a = a p h + a n a, 
where a ph = (M/T c )(dT c /dVt)(dVt/dM) is the usual (BCS) phonon contribution (fi is 
a characteristic phonon energy) and the non-adiabatic contribution is given by: 

n dT 

*na = J ¥ ^ , (34) 

where the parameter 7 = —M/n(dn/dM) has a weak logarithmic dependence on M 
(see Ref. ||). This parameter should not be confused with Euler's constant. Eq. (^) 
shows that the IC of T c depends on the doping of the conducting layer and on the 
relation T c {n). This result was used in Refs. |J[ Q to analyse the IE of high-temperature 
superconductors (see below). 



5.2 The Penetration Depth 

The concept of the non-adiabatic isotope effect relies on the fact that the charge- 
carrier concentration depends on the ionic mass because of the Jahn- Teller crossing 
of electronic terms. Obviously, any quantity that depends on n(M) will also display 
an unconventional isotope effect. One such quantity is the penetration depth of a 



magnetic field 5 given in Eq. flT7j). Note that the relation 5 2 ~ n s is also valid in the 
strong-coupling case (see, e.g., Refs. f|). From Eq. (H) one has 



1 5 dn s 8M 2n s dM [ } 

As mentioned in Sec. |2.7j , because of the relation n s = n(f(T), one has to distinguish 
two contributions to (3. There is a usual (BCS) phonon contribution, P ph , arising from 
the fact that (p(T/T c ) depends on ionic mass through the dependency of T c on the 
characteristic phonon frequency. This BCS contribution was discussed in Sec. |2|. In the 
present paper we focus on the non-trivial manifestation of isotopic substitution arising 
from the isotope dependence of the charge-carrier concentration n. Such an effect can 
even be observed in the low-temperature region where <p ~ 1. From Eq. (|35D and the 
relation n s = mp(T) it follows that 

R-R UR M & ^ ■ M 9n 

P Pph + Pna m + , W 



where n(M) is the normal-state charge-carrier concentration. Comparing Eq. fl3l|) and 
the second term of the right hand side of Eq. ( |36l) one infers that f3 na = —7/2 and thus 
establish a relation between the non-adiabatic isotope coefficients of T c and 5: 

n dT c 

This result holds for London superconductors. The equation contains only measurable 
quantities and can thus be verified experimentally. It is interesting to note that (3 na 
and ot na have opposite signs when dT c /dn > (which corresponds to the underdoped 
region of high-T c materials) . 

Note, finally, that in the presence of magnetic impurities the relation n s = nip(T) 
remains valid, but f(T) depends now on the direct scattering amplitude T 2 defined in 
Sec. |3]3| (this results, e.g., in the inequality n s {T = 0) < n in the gapless regime). As 
a consequence, magnetic impurities affect the first term in Eq. (^) (that now depends 
on T 2 ), but leaves (3 na and thus Eq. (|37D unchanged. 



6 OXYGEN ISOTOPE EFFECT 
IN HIGH-T C MATERIALS 

In this section we briefly discuss experiments done on high-temperature superconduc- 
tors in the light of the theory exposed in sections |3] and ^. A detailed study of the 
oxygen isotope effect in high-T c materials can be found in Refs. || f|, ^ f|. 

The oxygen isotope effect of T c has been observed in a number of experiments [|ZT|, 



72] , |73| . Here we discuss the results obtained on Pr-doped and oxygen-depleted 
YBa2Cu 3 7 _,5 (YPrBCO and YBCO respectively). We also discuss shortly the isotope 
effect of the penetration depth observed in La 2 _ x Sr :r Cu04 (LSCO). The analysis of the 
isotope effect in Zn-doped YBCO has already been presented in Sec. |3|. 

6.1 The Critical Temperature 

Let us first consider the isotope coefficient of T c for Pr-doped YBCO (YPrBCO). Several 
experiments have established that the Pr replaces Y which is located between the two 



CuC>2 planes of a unit cell. The doping affects YBCO mainly in two ways. First, 
it was shown that because of the mixed-valence state of Pr, holes are depleted from 
the CuC>2 planes (see, e.g., Ref. WW). Secondly, as for Zn substitution, praesodimium 



changes the magnetic impurity concentration of the system f56|, [F4j| . Regarding the IE, 
these experimental facts imply that through the first effect, the non-adiabatic channel is 
activated, whereas the second effect leads us to consider the magnetic impurity channel 
for the calculation of a. The IE is therefore described by a t = a p h + a m+na , where 
a p h = 0.025 |57j is the phonon contribution, and a m+na is defined by Eq. ([Tj^) with 
a given by Eq. fl3~4]) . The resulting expression depends on two parameters 7 and T s , 
characterizing the non-adiabatic and magnetic impurity channels, respectively. The 
latter quantity has been extracted from a fit to the relation T c (x) (x describes the 
Pr doping) J75| and is T s = 123K. The parameter 7 = 0.16 is determined from the 
mean-square fit to the experimental data. The result is shown in Fig. 8. The theory 
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Figure 8. Dependence of the isotope coefficient a on T c for Yi_ x Pr a; Ba2Cu307_ ( 5. Theory: 
solid line with 7 = 0.16, T s = 123K, oo = 0.025; Experiment: dc magnetization, resistivity 
and ac susceptibility from Refs. |7l], [56]]. 

is in good agreement with the experimental data. Note that the negative curvature at 
high T c 's reflects the influence of the non-adiabatic channel (magnetic impurities give 
a contribution with opposite curvature), whereas the positive curvature seen at low T c 
is due to the presence of magnetic impurities||. 

Oxygen-depleted YBCO is very similar to the previous case in that oxygen depletion 
both introduces magnetic impurities into the system and removes holes from the CUO2 
planes. The same equations as above can thus be used, however, with appropriate 
values of the parameters. The result is shown in Fig. 9, together with the experimental 
data from Ref. f57|. Note the sharp drop observed near T c = Q0K. This drop and the 
peak above it are related to the presence of a plateau in the dependence T c (n) near 
T, 



60K (see, e.g., Ref. |57|). Since Eq. (|J) contains the derivative dT c /dn, a na will 



be nearly zero at 60K. As one increases T c , the slope of T c (n) jumps to a large value, and 
decreases again to zero as T c — > 90K. This behavior of T c (n) explains the maximum 
appearing in a na (T c ) above 60K. The peak structure is thus specific to oxygen-depleted 
YBCO because neither YPrBCO nor YBCZnO display such a plateau in the relation 
T c (n). It would be interesting to verify this result, by measuring the oxygen isotope 
effect between 60 and 90K. 




The main difference between Pr-substituted (Fig. 8) and oxygen-depleted YBCO 
(Fig. 9) lies in the fact that the doping affects different ions of the system. Magnetic im- 
purities are introduced at different sites (on Y for the first material and in the chains or 
the apical oxygen position for the second system). Furthermore, non-adiabatic charge- 
transfer processes involve mainly Pr as well as Cu and O of the planes in YPrBCO, 
whereas they involve mainly the apical O in oxygen-depleted YBCO (to a lower extent, 
chain and in plane O and Cu are also involved). 

Contrary to YPrBCO, it is not possible to extract the value of T s (the contribution 
of magnetic impurities). The best fit to the few data gives 7 = 0.28 and T s = 15K. 
Two important points have to be noted concerning these values. First, because of the 
limited data available, the values may vary, although the order of magnitude will remain 
(see also Ref.||). The other remark concerns the value of T s for YPrBCO and YBCO. 
From these values one concludes that oxygen depletion introduces a smaller amount of 
magnetic moments than Pr doping. It would be interesting to perform more isotope 
effect experiments on oxygen-depleted YBCO and to determine the effective magnetic 
moment per depleted oxygen through other means. 



6.2 The Penetration Depth 

The only experimental observation of the isotopic shift of the penetration depth has 
been done on La2- x Sr x Cu04 (LSCO)Jft|. There are no data available on YBCO-related 
materials. Since only Pr-doped YBCO has one free parameter (7) we present the results 
of our theoretical calculations only in this case. The oxygen-depleted case was studied 
as a function of the two parameters 7 and T s in Ref. M. 

Let us begin with the case of LSCO. This case is simpler to study than YBCO 
related materials, since no significant amount of magnetic impurities has been detected 
in this material. Our analysis can thus be carried out with Eqs. (|34|) or ( j37|) . The 
isotope coefficient has been measured for Sr concentrations near x ~ 0.11 and x ~ 
0.15. The first concentration corresponds to the region where T c (x) experiences a small 
dip. The origin of this dip is not well-established but is probably related to electronic 
inhomogeneities and structural instability. Since many factors affect T c at this Sr 



concentration, it is difficult to interpret correctly the isotope effect of T c and 5. 

The second concentration at which an isotope shift of 5 has been measured is at 
optimal doping (T c is maximal). The experimental shift A5/5 ~ 2% and Eq. ([3]) allows 
one to determine the value of [3 na ~ 0.16 and thus 7 ~ 0.32 through Eq. fl37l) . This 
result is in good agreement with the values obtained for Pr-doped and O-depleted 
YBCO (see Ref. [[| for a discussion on LSCO). 

Let us now turn to the isotope coefficient of the penetration depth for YBCO-related 



materials. From the evaluation of the parameter 7 (see Eq. fl34|) and above) and from 
Eq. (|37|) one obtains the non-adiabatic contribution to the IC of the penetration depth 
in Pr-doped and O-depleted YBCO. One obtains (3 na = —0.08 in the first case and 
Pna = —0.14 in the second. 

Using Eqs. (0), (0) and (|34"D, as well as the values of the parameters 7 and T s 
given above, one can calculate the IC resulting from the contributions of magnetic 
impurities and non-adiabaticity for YPrBCO near T c . The result is shown in Fig. 10 
for three different temperatures. One notes that contrary to the case of magnetic 
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Figure 10. Dependence of the isotope coefficient j3 m+na on T c for T/T c = 0.75 (solid 

line), 0.85 (dashed), 0.95 (dotted). a = 0.025, 7 = 0.16, and f s = 123K (parameters 
for Yi_ a .Pr :c Ba2Cu307_ ( 5). 



impurities alone (one sets 7 = 0) studied in Sec. [I] (Figs. 2 and 3), the isotope effect 
of the penetration depth does not have the same qualitative behaviour as the IC of T c 
(compare Figs. 9 and 10). The change of curvature does not take place at the same 
value of T c (in the case of (3 m+na the change occurs at high T c 's and is barely visible 
on the figure; see Ref. || for a discussion of this point). One should measure the IC 
of 6 for YPrBCO and YBCO, since it would allow one to give a better estimate of the 
parameters and to test the theory. 



7 CONCLUSIONS 

We have reviewed different aspects of the isotope effect. In a first part, we have sum- 
marized the main theories developped to explain the isotope effect in conventional 
superconductors and we have discussed their relevance for high-T c oxides. In particu- 
lar, we have considered the effect of the Coulomb interaction, band structure (van Hove 



singularities), multiatomic compounds, anharmonicity and non-phononic mechanisms 
on the value of the isotope coefficient. These effects can account for most experimental 
data obtained on conventional superconductors. For example, they allow one to de- 
scribe the deviation from the standard BCS value a = 0.5 in transition metals or the 
inverse isotope effect (a < 0) found in PdH. The results show clearly that a vanishing 
or even a negative isotope coefficient can be obtained within conventional supercon- 
ductivity where the pairing between charge-carriers is mediated by the electron-phonon 
interaction. Generally, however, the small isotope coefficient is obtained only in con- 
junction with low T c 's. On the other hand, one can also obtain values of the isotope 
coefficient that are larger than a = 0.5 when anharmonic or band-structure effects are 
present in the system. 

In a second part, we have analyzed the effect of magnetic impurities, proximity 
contacts and non-adiabatic charge-transfer processes on the isotope coefficient (see also 
Refs. 0, f|, |^, H |7|). An important common feature of these factors is that they are 
not related to the pairing mechanism but, nevertheless, strongly affect the isotope 
effect. On notes that a > 0.5 is allowed in all three channels. Furthermore, these 
factors can induce a non-trivial isotopic shift of quantities such as the penetration 
depth 5. In the case of magnetic impurities and the proximity effect, the isotope 
effect of the penetration depth is temperature dependent. This phenomenon has not 
been investigated experimentally yet. In the presence of non-adiabaticity, we have 
established a relation between the isotope shift of T c and 5 for London superconductors. 

The theory presented in the second part of this review allowed us to describe the 
unconventional behavior of the isotope coefficient in YBa2Cus07 related systems (Zn 
and Pr-substituted as well as oxygen-depleted materials). The case of Zn-substituted 
YBCO was described by involving solely the magnetic impurity channel and did not 
require any fitting parameter. The theoretical curve also applies to conventional super- 
conductors. 

Our calculations suggest several experiments both on conventional and high- 
temperature superconductors. In particular, it would be interesting to measure the 
change of the isotope coefficient induced by a proximity system or magnetic impuri- 
ties in conventional superconductors to test our theory. Furthermore, more experiments 
have to be done on high-temperature superconductors so as to determine more precisely 
the parameters appearing in the theory. 
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9 APPENDIX 



We derive the explicit form of the function Rq given in Eq. ( pip (see Ref. M). The 
penetration depth calculated by Skalski et al. at zero temperature is given by 5~ 2 = 
-(Anne 2 /mc 2 )k(uj = 0, q = 0) with|0] 



^(0,0) 



I + Tfj- 3 



7T 



m 

R(fj) 



~-3 



3" 



71 



(38) 



for r, 77 < 1 (with f(rj) = arcosr,) or T < 1, 77 > 1 (with /(t?) = arcoshr,) and 
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(39) 
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7T 
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for r,77 > 1. We have introduced the notation T = T s /A, 77 — rjT — T 2 /A, 77 = T 2 /T s 
R(x) = 



-x 2 \ withx = T,fj and ft = [(T - l)(fj - 1)/(T + l)(fj+ 1)] 1/2 . A = A(T = 
0, r a ) is the order parameter in the presence of magnetic impurities. Eqs. (|38|) , fl39^) 
are valid when T s <C T 2 . These two scattering amplitudes, T 2 and T s (T s = T\ — T 2 ), 
defined by Abrikosov and Gor'kov|[49||, describe the direct and exchange scattering, 
respectively. One can calculate the magnetic impurity contribution to the IC at T = 
from Eq. (D) in a straightforward way using Eqs. ( 3§D and (|3~9"1). The result can be 
written as 



Pm(T = 0) 



«a Ki + K 2 



2 K(0,0) 

where a a is defined below, K(0, 0) is given by Eqs. (|38l), 
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for r, 77 < 1 (upper sign, f(rj) = arcos??) or F < 1, 77 > 1 (lower sign, /(t?) 
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for r, 77 > 1 and AT C /T C (T - 1) < 1. The last condition expresses the fact that the 
calculation is not valid in the immediate vicinity of T — 1. Eq. ( |4"0"D contains a a which 
is the IC of the order parameter A. In strong-coupling systems, «a has to be calculated 
numerically using Eliashberg's equations. Here we calculate the IC in the framework 
of the BCS model where «a can be calculated analytically. Indeed, from the relations 
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derived by Abrikosov and Gor'kov [A = A(T = 0, T s = 0) is the order parameter in 
the absence of magnetic impurities] one obtains 
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2 2f 



r < i 
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(44) 



In the BCS approximation one further has a& = ocq, where the last quantity was 
defined before as the IC of T c0 , that is, in the absence of magnetic impurities. 
Finally, one obtains Eq. (EBT) with R given by: 



R 



a a K\ + K 2 
2a K(0,0) 



(45) 



and a a is given by by Eq. fl44|). 
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